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I present results of a next-to-leading order calculation of three jet 
production at hadron colliders. This calculation will have many appli- 
cations. In addition to computing three-jet observables (spectra, mass 
ft. distributions), this calculation permits the first next-to-leading order 

studies (at hadron colliders) of jet and event shape variables. 



1 Introduction 



X 
b. 

One of the difficulties in interpreting experimental results is in assessing 
the uncertainty to be associated with the theoretical calculation. In QED 
and the weak interactions, one generally has confidence in the accuracy of 
leading order (LO) calculations because the couplings are sufficiently weak 
that higher order corrections are small. In QCD, however, the coupling is 
quite strong and it is difficult to obtain a reliable estimate of the theoretical 
uncertainty. 

One typically characterizes theoretical uncertainty by the dependence on 
the renormalization scale \x. Since one doesn't actually know how to choose 
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/i or even a range of /x, the uncertainty associated with scale dependence is 
somewhat arbitrary. One motivation for performing next-to-leading order 
(NLO) calculations is to reduce the scale dependence associated with the 
calculation. 

However, this is not the only benefit of an NLO calculation. There are 
times when the LO calculation is a bad estimator of the physical process. It 
may be that leading order kinematics artificially forbids the most important 
physical process. It could also be that the NLO corrections are simply large. 
Even if the overall NLO correction is relatively small, there may be regions 
of phase space, where NLO corrections are large. It is only in those regions of 
phase space where the NLO corrections are well behaved (as determined by 
the ratio of the NLO to LO terms) that one has confidence in the reliability 
of the calculation and can begin to believe the uncertainty estimated from 
scale dependence and it is only when one has a reliable estimate of the 
theoretical uncertainty that comparisons to experiment are meaningful. 

2 Methods 

The NLO three jet calculation consists of two parts: two to three parton 
processes at one-loop (the virtual terms) and two to four parton processes 
(the real emission terms) at tree-level. Both of these contributions are in- 
frared singular; only the sum of the two is infrared finite and meaningful. 
The virtual contributions are infrared singular because of loop momenta go- 
ing on-shell. The real emission contributions are singular when two partons 
become collinear or when a gluon becomes very soft. The Kinoshita-Lee- 
Nauenberg theorem [|| guarantees that the infrared singularities cancel for 
sufficiently inclusive processes when the real and virtual contributions are 
combined. 

The parton sub-processes involved are gg — > ggg ||] , qq — > ggg |J , qq — > 
QQg 0| , and processes related to these by crossing symmetry, all computed 
to one-loop, and gg -> gggg, qq -> gggg, qq ->■ QQgg, and qq -> QQQ'Q' 
and the crossed processes computed at tree-level. 

In order to implement the kinematic cuts necessary to compare a calcu- 
lation to experimental data one must compute the cross section numerically. 
Thus, it is not sufficient to know that the singularities drop out in the end, 
we must find a way of canceling them before we start the calculation. Sev- 
eral different methods of implementing this infrared cancellation have been 
successfully employed in various NLO calculations. The method we use is 



the "subtraction improved" phase space slicing method |7j. Phase space slic- 
ing || [| uses a resolution criterion s m i n , which is a cut on the two parton 
invariant masses, 

s ij = 2E i E j (l- cos 9ij). (1) 

If partons i and j have sy > s m i n they are said to be resolved from one 
another. (Which is not to say that a jet clustering algorithm will not put 
them into the same jet.) If Sjj < s m i n partons i and j are said to be 
unresolvable. One advantage of the s m j n criterion is that it simultaneously 
regulates both soft (E^ — ► or Ej — > 0) and collinear (cos B^ — ► 1) emission. 
In the rearrangement of terms, the infrared region of phase space is where 
any two parton invariant mass is less than s m j n . These regions are sliced 
out of the full two-to-four body phase space, partially integrated and then 
added to the two-to-three body integral. 

Because the infrared integral is bounded by s m i n , both the two-to-three 
and two-to-four body integrations are logarithmically dependent on s m ; n . 
Since s m ; n is an arbitrary parameter the sum of the two contributions must 
be s m i n independent. Thus, we have rearranged the calculation, trading a 
cancellation of infrared poles for a cancellation of logarithms of s m j n . The 
demonstration of s m ; n independence implies that we have correctly imple- 
mented the infrared cancellation. 

3 Results 

The results shown below were computed for the following kinematic con- 
ditions: the pp center of mass energy was 1800 GeV; at least one jet was 
required to have more than 100 GeV of transverse energy {Ex)-, while two 
more jets were required to have more than 50 GeV of transverse energy. All 
three jets were required to lie in the pseudorapidity range —4.0 < r\j < 4.0. 

The first test of the calculation is to demonstrate s m ; n independence 
of the cross section. In figure [l], the next-to-leading order cross section is 
computed for sixteen values of s m i n between 1 and 40 GeV 2 . We see that 
the computed cross section is stable over a wide range of s m i n . The next- 
to-leading order calculations were all performed using CTEQ3M parton 
distributions. The renormalization and factorization scales were chosen to 
be /I = 100 GeV. 

For comparison, the leading order calculation (computed using the CTEQ3L 
parton distributions but all other parameters the same) is shown as a solid 
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Figure 1: Next-to-leading order three jet cross section vs. s m j n . The leading order 
cross section is shown as a solid line. 



line. We see that for this choice of parameters, the magnitude of the next- 
to-leading correction is small. 

I have also computed the transverse energy spectrum of the leading jet 
(.En). Figure ^ shows a first attempt to explore the scale dependence of the 
calculation. The kinematic cuts (and parton distributions) for the results 
shown in figure |2] are the same as in figure ffl, but the renormalization and 
factorization scales for the upper, middle and lower curves are chosen such 
that \ip = /j,r = En/2,En,2En respectively. For the next-to-leading 
order calculation, s m ; n was chosen to be ~ 20 GeV 2 . 

Comparing to the leading order results (shown as solid lines), we see 
that the next-to-leading corrections are indeed small, and that the scale 
dependence is substantially reduced. 

4 Applications 

The next-to-leading order calculation of three jet production will have a 
wide array of phenomenological applications. 

4.1 Measurement of a s 

It should be possible to extract a purely hadronic measurement of a s . One 
possibility for such a measurement would be a comparison of the three jet 
to two jet event rate. Since both processes are sensitive to all possible initial 
states at tree-level, a next-to-leading order comparison should be relatively 
free of bias from the parton distributions. Because the measurement will be 
simultaneously performed over a wide range of energy scales, the running 
of a s can be used to constrain the fits and enhance the precision of the 
combined measurement. 
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Figure 2: Transverse energy spectrum of the leading jet for different choices of 
scales. The next-to-leading order results are shown as points and the leading order 
results as solid lines. In each case, the upper points correspond to /i = -En/2, the 
middle points to [i = Eti and the lower points to fi — 2Eti- 

4.2 Study jet clustering algorithms 

Because there are up to four partons in the final state, as many as three 
partons can end up in a single jet. This makes the three jet calculation 
sensitive to the details of jet clustering algorithms. This sort of study in 
pure gluon production J/]] uncovered an infrared sensitivity in the commonly 
used iterative cone algorithms. 



4.3 Study jet structure and shape 

Because there can be three partons clustered into a single jet, this calculation 
will allow truly next-to-leading order studies of the energy distribution in 
jets. Studies of jet production in deep inelastic scattering [|| show that the 
next-to-leading order correction for this variable is substantial and agrees 
rather well with experimental measurements. 



4.4 Study event shape variables 

There has been a long history of studying event shape variables like Thrust 
at e + e~ colliders. These measurements challenge the ability of perturbative 
QCD to describe the data and provide a means (other than event rate) of 
obtaining a precise measurement of a s . It will be interesting to see if one 
can make a meaningful study of such variables at hadron colliders. 
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